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We investigate the vector and scalar form factors relevant for semileptonic decays using 
maximally twisted-mass fermions with two flavors of dynamical quarks (Nf — 2). The simu- 
lations cover pion masses as light as 260 MeV and four values of the lattice spacing, ranging 
from ~ 0.05 up to ~ 0.1 fm, which allow to compute directly, for the first time, the continuum 
limit for the vector form factor at zero-momentum transfer, /+(0). The preliminary result is 
/+ (0) = 0.9544 (68 S f fl f ), where the error is statistical only. We also extrapolate both form factors 
to the physical point and study their momentum dependence. Our results are in good agreement 
with those obtained from a dispersion analyses of the experimental data. Together with the form 
factors, we analyze the ratio of the leptonic decay constants fic/fx, by imposing the constraint 
coming from the Callan-Treiman theorem, obtaining at the physical point fa/fx = 1.190 (8 sM )- 
Combining our results for /+(0) and fx /fit with the experimental measurements of the leptonic 
and semilpetonic decay rates, and using the determination of \V m i\ from nuclear beta decays, 
we determine the values of the Cabibbo angle |V HJ | from both Kfy and Ka decays, obtaining 
\Vus\ Ka = 0.2266 (17) and \V m \ K o- = 0.2258 (16). 
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1. Introduction 

Weak hadron decays are very interesting processes because measuring the decay widths for 
such processes allows us to extract some of the fundamental parameters of the Standard Model, 
namely the entries of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [jl]]. The 
needed theoretical inputs from non-perturbative QCD are the form factors and decay constants 
which parametrize the hadronic matrix elements relevant for each decay. 

In the case of the K — > semileptonic decay, the matrix element of the weak vector current 
can be written in terms of two form factors, the vector, f+{q 2 ), and the scalar, fo(q 2 ), form factors: 

(K(p n ) |V \K(p K )) = (p K +p K - Af f + (q 2 ) + A M fo(q 2 ) , (1.1) 

where A' 1 = q^ (M\ — M 2 )/q 2 and q^ 1 = (pk — p%)^ is the 4-momentum transfer. This decay is 
relevant for the determination of the CKM matrix element |V„ iS | ~ sin$, where 6 is the Cabibbo 
angle. 

In this contribution we present a lattice study of the vector and the scalar form factors per- 
formed by using the gauge configurations generated by the European Twisted Mass Collaboration 
(ETMC) with Nf = 2 maximally twisted-mass fermions. We present our preliminary results com- 
ing from two different strategies. The first one is a direct update of the analysis of Ref. [||], which 
was based on simulations at two values of the lattice spacing (a ~ 0.086 and a ~ 0.068 fm). Here 
we employ the results of simulations performed at four values of the lattice spacing, ranging from 
a ~ 0.05 up to ~ 0. 1 fm. This allows us to compute for the first time, in a well controlled way, the 
continuum limit for the vector form factor at zero-momentum transfer, /+(0). The second strategy 
is a multi-combined fit of the q 2 , M % and a dependencies of both the vector and the scalar form 
factors. We also analyze, together with the form factors, the ratio of the leptonic decay constants 
Ik/ fit, by imposing the constraint coming from the Callan-Treiman (CT) theorem [Q]. In this way 
we determine the full momentum dependence of the form factors at the physical point, finding re- 
sults which agree nicely with those obtained from a recent dispersion analyses of the experimental 
data Our preliminary results for /+(0) and fx/ fit are 

/+(0) = 0.9544 (68**) , hi fit = 1-190 (8, fflf ) , (1.2) 

where the quoted errors are statistical only. A detailed analysis of the systematic uncertainties, 
including the estimate of the quenching effect of the strange quark, is in progress and final results 
will be presented in a forthcoming publication. 

Using the experimental information from Kg$ and K& decays [Q], together with the precise 
determination of \V u d\ from superallowed nuclear beta decays, \V u d\ = 0.97425 (22) we obtain 
the following determinations of the Cabibbo angle 

\V US \^ = 0.2266 (17), \V US \ K ® = 0.2258 (16), (1.3) 

in good agreement among each other. 

2. First strategy 

We have performed the calculations of all the relevant 2-point and 3-point conelation functions 
using the ETMC gauge configurations with Nf = 2 dynamical twisted-mass quarks [Q] generated 
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at four values of /3, namely the ensembles A 2 — A4 at j5 = 3.8 (a ~ 0.101 fm), B2 — Bj at /3 = 3.9 
(a ~ 0.086 fm), C 2 - C 3 at J3 = 4.05 (a ~ 0.068 fm) and D 2 at j8 = 4.2 (a ~ 0.054 fm). The pion 
mass M K ranges between ~ 260 MeV and ~ 580 MeV and the size L of our lattices guarantees 
that M n L is always larger than 4.0 except for the ensemble B 7 {M n L ~ 3.7). For each pion mass 
and lattice spacing we have used several values of the (bare) strange quark mass m s to allow for a 
smooth, local interpolation of our results to the physical value of m s (see Ref. [Q]). 

Following Ref. the momentum dependence of the calculated form factors is fitted using 
either a pole or a quadratic behavior in order to determine the values of /+ (0) at each simulated pion 
and kaon masses. Then our results for /+ (0) can be smoothly interpolated (by quadratic splines) in 
terms of the kaon mass at a reference value M% , obtained by fixing the combination (2Mg — M%) 
at its physical value. 

Discretization effects are found to scale linearly with a 2 , as shown in Fig. [j], in agreement with 
the expectation of (automatic) & (a) -improvement of the maximally twisted-mass formulation [||]. 

1.002 



0.998 - 



0.994 



0.990 



0.986 



— 1 — 1 — 1 — 1 — 1 — 1— 
-O-M ~ 400 MeV 
t> M ~ 480 MeV 




0.000 



0.004 



0.008 



0.012 



2 li 2 N 

a (fm ) 



Figure 1: Results for /+(0) versus a 2 forM n ~ 400 and 480 MeV at Mk = ' . The dashed and solid lines 
represent linear fits in a 2 . 



Thus, as in Ref. [0], we perform the chiral extrapolation starting from the SU(2) ChPT predic- 
tion at next-to-leading order (NLO) [[J and adding to it both a NNLO term proportional to in 
order to analyze all our data up to M % ~ 580 MeV, and a linear term in a in order to correct for 
lattice artifacts, namely 



/+(0) 



1-T 



3 Ml 



4 {Anff 



Mi 



f-^log -f +c + Ml + d + Ml + Da< 



(2.1) 



In Eq. ( J2. 1[ ) / is the pion decay constant in the SU(2) chiral limit, and F + , c + , d + are SU(2) low- 



energy constants (LEC's) functions of the strange quark mass m s l . The quality of the fit (2.1) 
applied to our data with M n < 580 MeV is shown in Fig. ||, where we also compare it with the 
continuum limit result of a pure NLO fit (i.e., with d + = 0) applied to our data with M K < 400 
MeV. 



The LEC c+ depends also on the renormalization scale fi, but the whole result ( 2. 1 ) is independent on fl. 
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Figure 2: Results for /+(0) versus at Mk = M r ^ for four values of the lattice spacing a. The dashed 



lines represent the SU(2 )fit (2.1) at each value of the lattice spacing, while the solid line is the same fit in the 



continuum limit. The dotted line is the fit (2.1) at NLO ( i.e., with d+ — 0) applied to our data with M n < 400 



MeV and evaluated in the continuum limit. The vertical line corresponds to M% hys = 135.0 MeV. 



After correcting for the effect of quenching the strange quark, as discussed in Ref. [^] (namely, 
/+ (0) — /+^(0) = —0.0058 (28)), our preliminary result for /+ (0) at the physical point is 



/+(0) =0.9542 (6(W) 



(2.2) 



where the quoted error is statistical only. A final result including the estimate of the systematic 
uncertainties due to residual discretization effects, momentum and quark mass extrapolation, and 
the quenching effect of the strange quark, will be presented in a forthcoming publication. The 
preliminary value obtained for /+(0) agrees very well with our previous result of Ref. as well 
as with both the Leutwyler-Roos result flio| ] and previous lattice calculations performed with N/ = 
[0, N f = 2 ^ ||, and N f = 2 + 1 ^ dynamical flavours. 



3. Second strategy 

The second analysis consists of performing a multi-combined fit of the q 2 ,M % and a dependen- 
cies of the form factors in order to predict, at the physical point, both the vector and the scalar form 
factors not only at q 2 = 0, but also in the entire q 2 -region spanned by the experiments, i.e. from 
q 1 = to q 2 = q 2 max = {M K - M n ) 2 . 

We consider, for the vector and the scalar form factors, the following functional forms 

f2 



f+(s) = F + (s){\+C + (s)x + 



f (s) = F Q ( S )\\+C (s)x + 



Mi 



(4nf) 2 
M 2 



(4nf) 2 



- -xlogx — xTj (s) — Tj (s) 
lx\ogx+xT?( S )-T°(s) 



(3.1) 



where x = M 2 /M\ and s = q 2 jM\. The terms in the square brackets in Eq. (3.1) are derived from 
the NLO SU(3) ChPT predictions [ITS] for the kaon and pion loop contributions to the form factors 
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expanded in powers of x, keeping only the 
2 + (s) are then given by 



(jclogx) and ^(log(l — 5 , ))terms. The functions 



T, + (s) = [(l-s)\og(l-s) + s(l-s/2)]3(l+s)/4s 2 , 

T+(s) = [(\-s)log(\-s) + S (l-s/2)](\- S ) 2 /4s 2 , 

Z?(,y) = [log(l-^)+ 5 (l+V2)](9 + 7s 2 )/4/, 

Z2O) = [Cl-s)log(\-s) + s(l-s/2)}(\-s)(3 + 5s)/4s 2 . 



(3.2) 



It can be seen that the coefficients of the pion chiral log in Eq. (|3JJ) are in agreement with those 
predicted by SU(2) ChPT, both at q 2 = and q 2 = q 2 nax [§]. At q 2 = the leading chiral log has 
the coefficient (—3/4), while close to q 2 = q 2 max , i.e. for s ~ (1 — \/x) 2 , the functions 7\° 2 (j) also 
contribute to the chiral log leading, for fo(s), to an overall coefficient equal to (—11/4). 



The functions Fq :+ (s) and Cb,+ (s) in Eq. ( |3.l| ) are not predicted by SU(2) chiral symmetry. 
For Fq(s), we include in our analysis the constraint coming from the CT theorem which states 
that the scalar form factor fo{q 2 ) at the (unphysical) CT point, defined as q 2 CT = M\ —M 2 , differs 
from the ratio of the leptonic decay constants fy/fn by terms which are proportional to the light 
quark masses, namely: fo(q 2 = M^ — M 2 ) = f K / ' f n + &(m u ^)- Therefore, in the SU(2) chiral limit 



the scalar form factor fo(q ) at q CT 



Qmax 



M\ coincides with the ratio of the leptonic decay 



constants. Since the SU(2) chiral expansion of / f n is given at NLO by 



f_K = ft 
fn f 



l+Bx + 



Mi 



(4k f) 2 4 



-xlogx 



(3.3) 



where fg is the SU(2) chiral limit of f K , the CT theorem is equivalent to impose on F (s) the 
constraint 



F {s = \)=f K /f . 



(3.4) 



Inspired by the vector-meson dominance, we then adopt a pole behavior for Fq_ + (s) and a polyno- 
mial (quadratic) behavior for Cq,+ (s), namely 



Fo,+(s) =F/[\ -Ao,+j], C + (s) 



■ C ~\~ Cq _|_ s ~\~ Cq _j_ s 



(3.5) 



Finally, we take into account discretization effects by adding linear terms in a 2 to the parameter 



F and to the slopes Ao,+ of Eq. (3.5), as well as in the chiral expansion (3_3) of fa/ fn- As in the 
case of the first strategy, when we include in the fit all our data up to M n ~ 580 MeV, we also add 
in Eq. ([O]) a NNLO term of the form D 0)+ (s) x 2 , by expanding D 0j+ (s) = D + D^]_s + D^\s 2 . 

Our analysis involves a total of 120 data points with 16 free parameters, and we obtain a good 
quality fit with % 2 jd.o.f ' . ~ 0.8. The momentum dependence of the vector and scalar form factors 
extrapolated (for the first time) at the physical point is shown in Fig. |5[ The lattice results are also 
compared in the plots with those obtained from a dispersive fit of the experimental data [Q] from 
KLOE, KTeV, NA48 (without muons branching ratios) and ISTRA+, based on the parametrization 



of Ref. [17]. It can be clearly seen that our results are in good agreement with the data, in the whole 
range of q 2 spanned by the experiments. 
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Figure 3: The q -dependence of the scalar (a) and vector (b) form factors. The dashed lines correspond to 



the region allowed at 1(7 by the fit of the lattice data according to Eqs. (3.1) and (3.3). The shaded blue area 
is the 1(7 region obtained from a dispersive fit /|7^/ of the KLOE, KTeV, NA48 (without muons branching 
ratios) and ISTRA + data performed in Ref. /Q/ . 



The preliminary results of this second strategy for the form factor at zero momentum transfer 
/+(0) and the ratio f K /f n are 

/+(0) = 0.9546 (76. VM ) , f K /f n = 1.190 (8, fflf ) , (3.6) 

where, as in the case of the first strategy, the central value of / + (0) obtained from the fit has been 
shifted by 8f + = /+(0) — f+®(0) = —0.0058 in order to correct for the effect of quenching the 
strange quark. 



4. Results and conclusions 



Our preliminary best result for /+(0) is obtained by averaging the determinations (2.2) and 



(3.6) of the form factors obtained from the first and the second strategies respectively, leading to: 

/ + (0) =0.9544 (68„ flf ). (4.1) 



Eqs ( P^q ) also provides our estimate of the ratio fK/fn obtained by fitting the lattice data for the 
decay constants together with the semileptonic form factors and imposing the contraint coming 
from the CT theorem: 

f K /fn = 1.190(8**). (4.2) 

A careful analysis of the systematic uncertainties in the present calculation is still in progress. The 
results in Eqs. Q4.1[ ) and ( |4.2[ ) are in good agreement with the previous ETMC determination of the 
form factor, /+(0) = 0.9560(84) 0, based on simulations at only two values of the lattice spacing, 



and with the more extensive analysis of the meson decay constants presented in Ref [18], which 
quoted/*//, = 1.210(18). 



6 



Improved analysis of the scalar and vector form factors ofkaon semileptonic decays 



L. Orifici 



Combining our results ( |Q| ) and ( |4.2[ ) with the latest experimental averages [Q] \V US \ ■ /+(0) = 
0.2163 (5) from Ka decays and |V„. S /V M ^| • fx/fx = 0.2758 (5) from decays, and with the 
determination \V u d\ = 0.97425 (22) from nuclear beta decays [§], we obtain for the Cabibbo angle 
the values 



V us \ Ka = 0.2266 (17) , \V us \ Kn = 0.2258 (16) . (4.3) 



The average of the and Ka results in Eq. (4.3) can be further combined with V ll( ] 



0.97425 (22) and V ub = 0.00376 (20) Q19p to test the unitarity of the first row of the CKM ma- 
trix, for which we find 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1.0003 (8) . (4.4) 
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